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Abstract. The renormalized photon and electron propagators are expanded over planar binary trees. Ex- 
plicit recurrence solutions are given for the terms of these expansions. In the case of massless Quantum 
Electrodynamics (QED), the relation between renormalized and bare expansions is given in terms of a Hopf 
algebra structure. For massive quenched QED, the relation between renormalized and bare expansions is 
given explicitly. 

PACS. 12.20.-m Quantum electrodynamics - ll.10.Gh Renormalization 



1 Introduction 

Wightman commented jjj : "Renormalization Theory has a 
history of egregious errors by distinguished savants. It has 
a justified reputation of perversity; a method that works 
up to 13-th order in the perturbation series fails in the 14- 
th order." Although renormalization theory is often con- 
sidered to be well understood, it is still a difficult subject 
plagued with considerable combinatorial complexity. 

However, renormalization is not just a recipe to extract 
a finite part from an infinite integral, since it has a deep 
physical meaning. It was a guide to elaborate the theories 
of weak and strong interactions. It can be used to build 
consistent Lagrangians: the minimal coupling Lagrangian 
of scalar electrodynamics misses a quartic term which is 
reintroduced by renormalization Q. Renormalization is 
also linked to the irreversibility of the macroscopic Uni- 
verse [|| . Other arguments in favor of renormalization have 
been given by Jackiw Moreover, thanks to the work 
of Kreimer and Connes ||,[|0,[8,|^], renormalization has 
become a beautiful mathematical theory. 

There is a class of physicists who think that a quantum 
field theory is entirely contained in its Feynman diagrams. 
Since two members of this class, Veltman and 't Hooft, 
have been recently awarded the Nobel prize in physics, it 
seems rather out of place to try to formulate quantum field 
theory without Feynman diagrams. Nevertheless, we shall 
pursue our attempt to build such a theory by studying how 
renormalization can be implemented in the framework of 
planar binary trees. 

The typical equations of QED are 'implicit' functional 
derivative equations. In |nj, the tree expansion method 
enabled us to establish an explicit recursive solution for 
the photon and electron propagators. For renormalized 



QED including mass renormalization the functional equa- 
tions change drastically, and it is not a priori clear that 
a similar explicit recursive solution can be given. In this 
paper, we obtain an explicit recursive solution for renor- 
malized QED. In other words, the terms of the tree ex- 
pansion for the renormalized photon and electron prop- 
agators are written as an integral of renormalized terms 
for smaller trees. In the case of quenched QED, this re- 
cursive equation is transformed into an explicit relation 
between the renormalized and bare terms. Finally, we ex- 
hibit a Hopf algebra that encodes the renormalization of 
massless QED. This Hopf algebra is neither commutative 
nor cocommutative. 

The plan of the paper is the following. We first intro- 
duce the problem of renormalization from different points 
of view, then we give the essential equations for the renor- 
malization of QED. The tree expansion method is pre- 
sented in detail, and used to derive the recursive equa- 
tions for the renormalized electron and photon propaga- 
tors. Then detailed relations between the renormalized 
and bare propagators of quenched QED are given. Finally 
the Hopf algebra of massless QED is described. A first ap- 
pendix gives some proofs, a second one gives the relation 
between renormalized and bare propagators for massless 
QED. A last appendix expresses the smallest planar bi- 
nary trees as a sum of Feynman diagrams. 

A planar binary tree can be considered as a sum of 
Feynman diagrams (see appendix 3). To some extent, this 
sum is natural. For instance, all the diagrams obtained by 
permutations of the photon external lines of a Feynman di- 
agrams cannot be distinguished by an experiment. There- 
fore, it is natural to regroup all these diagrams under a 
single symbol. However, when renormalization comes into 
play, each Feynman diagrams is renormalized by a num- 
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ber of Feynman subdiagrams multiplied by scalar coun- 
terterms. It is not a priori clear that these subdiagrams 
and counterterms can be grouped into the same trees as 
the original Feynman diagrams. The main point of this 
paper is to show that this is possible for QED. In other 
words, we provide an algebraic structure on planar binary 
trees which is compatible with renormalization. 

In this paper, we consider only the renormalization 
of ultraviolet divergences, and we assume that infrared 
divergences are regularized, for instance by introducing a 
photon mass. 



2 Renormalization 

An enormous literature has been devoted to the renormal- 
ization theory. The reader is deferred to [[11,12,1^] and 
references therein for historical and conceptual aspects of 
renormalization. Here we shall concentrate on its technical 
aspects. Renormalization theory can be considered from 
at least three different points of view: the Dyson-Salam 
method, the extension of distributions and the product of 
distributions. 



2.1 The Dyson point of view 

According to the first point of view, perturbative quantum 
field theory yields divergent integrals in Fourier space, and 
renormalization is a technique intended to extract a finite 
part from them. A picture of how this could be achieved 
was first given by Dyson in 1949 @ and Salam fT|,[l6| 
in 1951. Explicit formulas were proposed by Bogoliubov 
and Parasiuk [|l7j and finally proved by Hepp pl||l9|. This 
method is general, in the sense that it can be used for any 
quantum field theory, whether renormalizable or not. 

To understand this renormalization process, it is very 
useful to treat a one-dimensional toy model of overlapping 
divergences proposed by Kreimer [0. Let 



f(x,y,c) = 



1 



y 



x+cx+yy+c 
We want to give a meaning to the integral 



/CO rC 
dx J 



dyf(x,y,c). 



Power counting is applied as follows. If we substitute 
Xx to x and take the limit A — > oo, we see that 1(c) varies 
as A = 1, and the integral is logarithmically divergent 
for x. Similarly, it is logarithmically divergent for y. If 
both x and y are multiplied by A, the integral 1(c) varies 
as A 1 in the limit A — > oo. Then, 1(c) is linearly divergent 
for the variables x, y. In the Dyson-Salam renormalization 
scheme, we first fix y in f(x,y,c), we keep the part of 
y, c) which does not depend on x (i.e. yj (y + c)) and 
we take the value of the rest (i.e. x/((x-\-c)(x-\-y))) at c = 
and y = (i.e. 1/x). The product of these two factors (i.e. 
y / (x(y J r c))) is called a counterterm and is subracted from 



f(x,y,c) to remove the logarithmic divergence for x. This 
procedure produces the term 



f(%,y,c) 



y 



y 



xy + xc + yc 



x(y -he) y + c x(x + c) (x + y) ' 



which is now convergent for the integral over x (it varies as 
A -1 by power counting). If we make the same subtraction 
while fixing the variable x we obtain, subtracting both 
counterterms: 



g(x,y,c) = f(x,y,c) 



y 



x(y + c) y(x + c) 



x 2 y 2 + xy 3 + yx 3 + cy 3 + cx 3 + cyx 2 + cxy 2 
x(x + c)(x + y)(y + c)y 

The result is disappointing, because g(x,y,c) is now lin- 
early divergent if x is multiplied by A, if y is multiplied by 
A and if x and y are both multiplied by A. In other words, 
g(x,y,c) is still more divergent than f(x,y,c). The mira- 
cle happens when we subtract the global linear divergence 
of g(x, y, c). The final term 



f(x, y, c) = g(x, y, c) - g(x, y, 0) - c 
2 xy + cx + cy 



dg(x,y,0) 
dc 



" x(x + c)(x + y)(y + c)y' 
is now absolutely convergent for x, for y and for x,y. 

2.2 The extension of distributions 

From a mathematical point of view, renormalization the- 
ory can be considered as the problem of extending a dis- 
tribution to a larger domain. 

The standard example is 1/x. If <j)(x) is a test function 
that vanishes at 0, then 



L 



exists. The question is how is it possible to extend this 
distribution to general test functions. The existence of this 
extension is ensured by the Hahn-Banach theorem J2(J and 
a formula for such extensions is 



/ 

J\x\<a 



dx 



4>{x) - 4,(0) 



J • 

J\x\>a 



(j)(x) 



dx 

l\x\>n X 



for any positive parameter a. Hence various extensions are 
possible, that are parametrized by a. Notice that the dif- 
ference between two such integrals for a and a f is (log a! — 
loga)0(O). Therefore, as distributions, two extensions of 
1/x differ by \og A5(x) for some A. The peculiarity of 
quantum field theory is that A can be determined by ex- 
periment. 

The mathematical conditions for the existence of such 
an extension were investigated by Malgrange [^lj and Estrada 
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This extension method can be used to calculate, in 
some cases, the product of two distributions. For instance, 
by Fourier transform, it can be shown that S(x — a)S(x) = 
(5(a), for a / 0. However, if a = 0, the Fourier transform 
of the product diverges. This is exactly the same type of 
divergence that is met in the usual presentation of renor- 
malization. The product of distributions S(x) 2 is zero for 
x 7^ 0, thus a possible extension is S(x) 2 = C5(x), where 
C is a constant determined by experiment. 

In quantum field theory, causality, Poincare invariance 
and unitarity were used by Stuekelberg and coll. to pro- 
vide a prescription to carry out this extension |23|,p4,p5|, 

tBogoliubov and coll. systematized this construction 
E^EpL^I , which took its final form with Epstein and 



Glaser [|30|. Nowadays, the extension method is called the 
"causal approach" , and the case of QED is treated in de- 
tail in Ref-lfS. 

A (correct) proof of the validity of Bogoliubov's method 
was finally given by Hepp in 1966 and by Zimmer- 
mann [Q in 1969. 

Recently, the causal approach has been reintepreted 
in terms of microlocal analysis |3^| . This enabled these 
authors to adapt the causal approach to quantum field 
theory in curved spacetime. 

An up to date and clear presentation of the causal 
approach can be found in Ref . Q . 

2.3 The product of distributions 

The most radical approach to renormalization would be 
to define a product of distributions, which could lead to 
a nonlinear theory of distributions. Schwartz has shown 
that this is impossible in general ]35| , but the notion of 
distribution can be extended to a more general kind of 
functions which can be multiplied. For a comparison with 
experimental results, we must project these new functions 
back onto standard distributions. 

This approach was investigated by various authors |36[ 
||,||,|3^g,|l|,||,|43[||. The main drawback of these new 
generalized functions is that they lead to very intricate 
calculations. For instance, it is not difficult to show that 



\ + tt 2 5(x) 2 , 



but the computation of (1/x) 3 is already intractable. To 
understand this striking identity, we start from the contin- 
uous function f(x) = x\og \x\ — x, which defines a distri- 
bution by J dxf(x)(j)(x) for any test function <j){x). Then 
the distribution 1/x is defined as cPf/dx 2 , the distribu- 
tion 1/x 2 is —d 3 f/dx 3 , and (1/x) 2 is the product of the 
distribution 1/x with itself. 

In spite of their complexity, these new generalized func- 
tions have found some applications in physics |I^,f47|,[48[ 
f49| . For instance, a definite value could be given to the 
curvature of a cone at its apex p7| . 

Notice that, as for the extension of distributions, mi- 
crolocal analysis is taking a growing importance in the 
study of the new generalized functions [M . 



3 Renormalization of QED 

QED was renormalized to all orders by Dyson We can 
now interpret his prescriptions in the framework of the 
Schwinger equations. It is standard to define free, bare 
and renormalized propagators. The free electron Green 
function S°(q) is the Green function for an electron with- 
out electromagnetic interaction. The bare electron Green 
function S(q) is the Green function for an electron with 
electromagetic interaction, but without renormalization. 
In the perturbation expansion of S(q), all terms (except 
the first one) are infinite. The renormalized Green function 
S(q) is the Green function for an electron with electro- 
magnetic interaction, after renormalization. Similarly we 
define D°(q), D(q) and D(q) as the free, bare and renor- 
malized photon Green functions. [] 



3.1 The free propagators 

The free electron propagator is 

S°(q) = (j-q-m + ie)- 1 . 
The scalar product is defined by 

A/j, 

where the pseudo-metric tensor is 



(1 
0-100 
0-10 

\0 -1, 



All the electron propagators S°(q), S(q) and S(q) are 
4x4 complex matrix functions of the 4- vector q. If / is 
the 2x2 identity matrix and a x , a y , <j z the Pauli matrices 



&x = 



1 

1 



(j v = 



-i 

1 



(7 Z = 



the Dirac matrices can be written 7 
a x \ 2 I (7 V 



1 
-1 

I 
-/ 



7 



-CF X 



7 



-dy 



cr z 
~(Jz 



The free photon propagator D°(q) is a complex 4x4 
matrix with components D^ u (q) defined by 



KM 



< 2 + ie 



+ (i-vo 



(q 2 + ie) 2 ' 



1 Strictly speaking, a propagator is a one-particle Green 
function, but in this paper propagator and Green func- 
tion will be used indifferently. For simplicity, fermions (elec- 
trons+positrons) are called electrons. 
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The term l/£ was introduced by Heisenberg to make D°(q) 
non singular. The Green function used in classical electro- 
dynamics is D 0T (q) defined as 



q 2 + ie ' (q 2 + ie) 



2 ' 



A tensor T^ v (q) such that q^T^ v (q) = is called trans- 
verse. It can be checked that D°^(q) is transverse, and non 
singular in the space of transverse tensors. 

Up to an eventual factor z, the expressions for the free 
Green functions S°(q) and D Q uv (q) are standard (see, e.g. 
Ref.g] p.93 and p.36, Ref.Q p. 184 and p. 190, Ref.[|| 
p. 218 and p. 253, for a complete description and a deriva- 
tion). 



3.2 The bare propagators 

The Schwinger equations for bare electron and photon 
propagators were given by Bogoliubov and Shirkov [ p9[ 
and transformed into the following integral equations in 



S(q) = S°(q) 



SS(q-p) 



(2?r)4 , — ^ eoSA O x/{p) 
D^(q)=D%(q) 



, (1) 



(2tt)4 l/ e (L4 A/ (-g) 



^av(^), (2) 



where (p) is an external electromagnetic field and j^u^ 

is the functional derivative evaluated at A^(p) = 0. 

The longitudinal part of the bare photon Green func- 
tion is not modified by the interaction g, and D\^(q) can 
be written as the sum of its transverse and its longitudinal 
parts: 



(3) 



where D% (q) is transverse. In Eq.(||), the photon propa- 
gator D\> u (q) is not integrated, it just multiplies the in- 
tegral. This is not very convenient and it will be useful to 
introduce the bare vacuum polarization, denoted II\^{q) 
and defined by 

\ D ~\M = ~ fdx^) - iq\q» + n X p(q). (4) 

The vacuum polarization tensor II\^(q) is transverse 
If we multiply (|4|) by (||), we obtain 

9x~^f= " Q 2 9X») + n^{q))D T ^{q). (5) 

The left-hand side of Eq.(||) is the projector onto the trans- 
verse tensors. 

We show in section |l3.2| that 



SS(p) 



e SA°(-q)- 



(6) 



3.3 Renormalized propagators 

To obtain the Schwinger equations for the renormalized 
propagators, the best is to start from the renormalized 
Lagrangian, and to follow the steps given by Bogoliubov 
and Shirkov p9j , Itzykson and Zuber Q or Rochev [ p3| . 
However, to give an idea of the result, we introduce some 
of Dyson's recipes. 

The longitudinal part of the photon Green function is 
not modified by renormalization , and the renormalized 
photon propagator can be decomposed as 



Dx,(q) = D^(q) - - 



qxQv 



(7) 



£ (q 2 +ie) 2 ' 

where D^(q) is transverse. 

Then, we introduce Dyson's relation between renor- 
malized and bare Green functions (Ref.(|] p. 414): 

S(q)Z 2 = 5(g), (8) 



Z 3 e 2 
m 



DUq), 



e , 
m 



5m, 



(9) 
(10) 
(11) 



where Z2 and Z3 are (infinite) scalars independent of 
and e is the renormalized charge. Equation (O) was con- 
jectured by Dyson [Q and proved by Ward (]2|, p. 413). 
Finally, the external field is renormalized as A\, so 
that 



e A° x = eA y 



(12) 



To introduce the mass renormalization, we must start 
from the differential form of the Schwinger equation for 
the bare electron propagator, where we reintroduce the 
external field, for later convenience, 



[ij-d-mo- e 7 ■ A°(x)]S(x, y; A ) = S(x - y) + 



iei / d"z^D w {x,z;A^) 



<5S(x,y;A°) 



Se A%z) 



■ (13) 



Dyson showed that relations (|8|-|12D, are valid in the 
presence of an external field. We use them in Eq.([L3|) to 
obtain 

[27 • d — m]S(x, y; A)Z2 = S(x — y) — 5mS(x, y\ A)Z2 



d*z^D w {x,z;A) 



5S(x,y;A) 



Zo. 



(14) 



SeA p (z) 

In Eq.(|l4|), we have changed the gauge parameter £o of 
D^ p (x, z; A ) into £ = Mo (see Ref.g p. 414). 
If we multiply Eq.(|l|) by 

S°(z, y;A) = [ij-d-m-ej- A}' 1 (15) 

and integrate over x we obtain the integral Schwinger 
equation for the renormalized electron propagator, 

S(x,y;A)Z 2 = S°(x,y;A) 



+ie 2 J d 4 zd 4 z'S°(x, z; A)j x D xx , (z, z'; A) 
-5m [ d i zS°(x,z;A)S(z,y;A)Z 2 . 



SS(z,y;A) 
e5A y {z>) ■ 



(16) 



Christian Brouder, Alessandra Frabetti: Renormalization of QED with planar binary trees 



5 



In Eq. ( |16| ) , we put A = and we Fourier transform to 
find 

S(q)Z 2 = S°(q) - SmS°(q)S(q)Z 2 

+ie s s>) /l^ ixvW |^ z , (17) 

This equations was given by Bogoliubov and Shirkov 
p9[ , as well as by Itzykson and Zuber ([[|, p. 481), ex- 
cept for the mass counterterm dm which was apparently 
overlooked by these authors. A complete derivation can 
be found in J53[ (notice that his Sm is our Z2S171). 

To obtain a convenient Schwinger equation for the 
renormalized photon propagator, we must introduce the 
renormalized vacuum polarization Ilx^q), defined by 

= " « 2 ^m) " ^ + n xM- (18) 

It may be useful to compare these definitions to those 
of Itzykson and Zuber ||: D^ v = -iG^, = -iu^. 
If we multiply (|l|) by (g), we obtain 

9x- S *f= ((^ ~ <rW) + n^{q))D T ^{q).{l9) 
If we compare Eqs.(|5|) and ([l9|), and use (^), we find 

+77 AM (g)). (20) 
Therefore, using Eqs. (||) and Jl^ ) 

Introducing this equation into (pp|), and using Eq.(|l0|) we 
obtain 

^(g) = (^3 - 1)(^A^ - ^W) 

+ ^ <2i) 

Equations (|l7| ) and ( ^T|) will be the bases of recur- 
sive expressions for the renormalized electron and photon 
propagators. 

4 Tree expansion of propagators 

For the convenience of the reader, and because the nota- 
tion of Ref. (HJ as been modified Q we recall the descrip- 
tion of photon and electron propagators in terms of planar 
binary trees. 

2 There are no longer trees with black or white roots. The 
color of the root is now indicated by the function ip itself. <p(t] q) 
corresponds to a tree with a black root, iffj,u(t; q) to a tree with 
a white root. This notation is more elegant than the one used 
in 0. 



4.1 Trees and propagators 

The main trick of Ref.[|l0| was to write each propagator 
as a sum indexed by planar binary trees, to be defined in 
the next section. 

The bare electron Green function in Fourier space, 
S(q) is written as a sum over planar binary trees t 

S(q)=J24 lt ^°(t;q). (22) 

t 

Here eo is the bare electron charge (i.e. the electron charge 
before renormalization). The fact that the expansion is 
over (and not eo) was justified in Ref.Jn]]. Similarly, 
the renormalized electron Green function is expanded over 
planar binary trees 

S(q) = Y,e 2ltl ?°(t;q). (23) 

t 

In Eq. (p3|) e is the renormalized (finite) electron charge. 

The bare and renormalized photon Green functions are 
written as 

t 

D, v {q)=Y, em ^U t ^)- ( 24 ) 

t 

For the renormalization of the photon Green func- 
tion and the vacuum polarization, it will be necessary 
to distinguish the photon Green function D jJjU (q) and the 
transverse photon Green function D^ u (q). Since all terms 
(p°^(t;q) and (f°^ u (t;q) are transverse for t ^ i, the trans- 
verse renormalized propagator is 

|t|>0 

where = D°£(q). 

The bare and renormalized vacuum polarization are 
expanded similarly: 

n x ,(q)=Y,4 lt ^U t ^)> ( 25 ) 

|t|>0 

nxM = E e 2|t| $ M (*;<z)- ( 26 ) 

|t|>0 

For later convenience, we finally define 
so that 

[D°~\(q) = (qxq» ~ q 2 9^) ~ ^ (27) 
= -&A«„. (28) 

A few identities will be useful in the sequel 

D%(q)(q^-q 2 g^ = -q 2 DZ(q), 
^Ar,q)D 0X '^(q)^^(r,q)=^%(r,q), 
D° xx ,(q)r X ' ti '(r,q)Dl, ll (q)=D x l(q). 
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Notice that —q 2 D x ^(q) is the projector onto the transverse 
tensors. 

The tree representation of the photon propagator en- 
joys the following property pj 

Vf P ;g) = V° A (- V^)[(Z>VT A Wav(^)- 

This equality is non trivial only if t\ i . But then 
(p ^ x (\\/t r ; q) and yP xly (ti\q) are transverse p)| , this cancels 
the term proportional to £ and we can write 

V%(tl V U- q) = (pl x ( , V t r - q)^ XX ' ( , ; q)^ u (t i; q) (29) 

The fact that tp\ v (t\ V t r ;q) is a product for t\ ^ i can be 
checked in appendix 3. 



As shown in section 13.2, ip® u (ti Vt r ) = if ti ^ i and 
i>%{ >Vt r ;q) = -V£ A ( i ; <?V° AA ' ( i V t P ; g)^v( ' ; «)(30) 



4.2 Trees and renormalization constants 

According to Dyson's multiplicative renormalization of 
QED, we define three renormalization constants Z2, Z3 
and 5m. We expand these constants over planar binary 
trees 

Z 2 = Y / e 2l %(t), with C 2 (0 = 1, (31) 
t 

Z 3 = ^e 2 \%(t), with Cs(i) = 1, (32) 

t 

5m = ^2e 2 ^C m (t), with Cm(0 = 0. (33) 
t 

These expansions will be the basis of the tree by tree 
renormalization of QED. 

4.3 Planar binary trees 

A planar binary tree is a tree with a designated vertex 
called the root. To follow the notation of Loday and Ronco 
J54j , we write the root vertex as 1 . The other vertices are 
not explicitly drawn, but they are at the ends of each edge, 
which are 

or /. The trees are binary because each vertex has either 
zero or two children. They are planar because ^ is differ- 
ent from f. The planar binary trees have an odd number 
of vertices and for each tree t we define |t| as the integer 
such that t has 2\t\ + 1 vertices. In other words, \t\ is the 
number of internal vertices. We call Y n the set of planar 
binary trees t with such that \t\ = n. 

The planar binary trees with up to 7 vertices are 

Yo = {>}, 

Yi = m, 

y 2 = o,?}, 

y 3 = {>^,Y^f } . 



We denote Y the set of all planar binary trees 

00 

Y= \jY n . 



n=0 



Finally we consider the operation of grafting two trees, 
V : Y p x Y q — ► Yp + q + i, by which the roots of two trees t\ 
and £2 are joined into a new vertex that becomes the root 
of the tree t = t\ V £2, cf IM- For instance 



Y V Y = Y. 



(34) 



It is clear that any tree t, except the 0-tree 1 , is the grafting 
of two uniquely determined trees t\ and t r with orders 

\u\Au\ < \t\-i. 



4.4 Recursive equations for bare propagators 

In Ref . pi , we have obtained recursive relations for ip(t) 
and (fx^Jt). 

For the electron propagator, ip(t) satisfies the recursive 
relation 



k) 



^(t;q;{\,p} hn )=S (qh^ 
X(p n - 1 (t;q+p 1 ;{X,p} 2 , n ) 

+*E/ j^S (qh X ip k xx ,(tr,p;{X, P K 

Xip n E - k+l (t r ; q-p;\',p+ P k , {X,p} k+hn ), (35) 

where we have noted P k = Pi + ■ ■ ■ + Pk, (Po = 0) and 
{A,p}i,n = Xi,p\, . . . , X n ,p n . The initial data are 

V°(r,q) = S°(q), 
l p 1 (r,q;X 1 ,p 1 ) = S (qh^S°(q + Pl ), 

f n ; v, {a,p}i,„) = s°(«h Al s°(g + Pl ) 7 A2 • • • 7 A " 

S (q + P i + ---+Pn)- (36) 

The symbol <^ +1 (t; q; z, {z}i, n ) is defined as the sum 
of n terms, where the first variable z = (X,p) is exchanged 
in turn with all the variables Zi = (Xi,Pi). 

tp" +1 (t;q;z,{z} hn ) = tp n+1 (t;q;z,{z} hn ) 

+<f n+1 ( Zl ,z, {z} 2 , n ) + ■■■ + <^™ +1 ({^}i, n _i, z, z n ) 
+p n+1 ({z} hn ,z). 

For the photon propagator, ^p^ v (t) satisfies the recur- 
sive relation 

n p n4 

<P%{t;q;{\,p}i, n ) = E / (2^^ 

xtr[ 7 A ^ +1 (rr;p; A', -q - P k , {X,p} hk )] 
x V n > 7 v k {t l -q + P k -{X,p} k+hn ), (37) 
with the initial data 

<pl v (r,q) = DP v (q), 
^(i;«;{A,p}i, n ) =0 for n>l. 
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In this paper, we use a non symmetrized (f n (\). How- 
ever, for the validity of Ward-Takahashi identities, it would 
be necessary to handle symmetrized expressions 

— ^2 ^ n ('5^5 . . . , A a ( n ),_p a ( n )). 

Such expressions seem complicated. It is possible that re- 
cent string-like methods help handling them pq,p6|. 



4.5 The pruning operator 

In this section, we introduce the pruning operator P which 
will prove very useful to obtain a recursive expression for 
renormalized propagators. If t is a tree, P(t) is a sum of 
n(t) terms of the form Uj <S>Vj, where Uj and Vj are planar 
binary trees. More mathematically 



'3 ^ u 3 • 



(38) 



i=i 



Before we fully define P(t), we want to show why it is 
useful. If, for each tree t, cp(t) and ip(t) are 4x4 complex 
matrices, we call the convolution of cp and i(j the quantity 

n{t) 

(^★^)(t) = ^2<p(Ui)l/)(Vi). 

2=1 

The main property of this convolution was established in 
Ref.0. If 

X(A) = £ \W x (t) and Y(X) = ^ \ w y(t), 
t t 

with x(\) = y(\) = 0, then 

x(\)Y(\) = ^2^ tl (x*y)(t). 



In other words, the pruning operator and the convolution 
enable us to multiply series indexed by planar binary trees. 

This nice property justifies the trouble of introducing 
P(t). First, n(i), the number of terms in Eq.(^8|), is denned 
by n( i ) =0 and 

n(t) =0 if t = ti V i , 

n(t) = n(t r ) + 1 if t = tiVt r , t r +\. (39) 
Finally, P(t) is determined recursively by P( i ) = and 
P(t) =0 if £ = ti V i , 

n(t r ) 

= (t« v i ) © t r + ( f i v ® 

if t = ttVt r , t r ^ i. (40) 



The trees and i>j in Eq. (f40|) are generated by Eq.(38) 
for t = t r . For instance 

p(,) = p(y) = P( N ?) =P(>) = P(^r) = o, 

P(f) =Y©Y, P(Y)= ^ ©Y, 

p(^) = y© N t, p(f) =T®r+r®r. 

We show in the appendix that the pruning operator is 
coassociative, that is 



(P ®id) © P = (id P) P. 



(41) 



Therefore the convolution is associative. 

We consider on trees the structure of an associative al- 
gebra T(Y) given by the (non commutative) tensor prod- 
uct, T(Y) = Y © Y m © Y® 3 © . . . , and we set the root i 
as the unit: i © t = t © i = t. Then we extend P to T(y) 
as a multiplicative map, but P does not preserve the unit, 
since P( i ) is not equal to i © i . We can define a coproduct 



> t + P{t) + t < 



This zA p is the coproduct of a Hopf algebra over planar 
binary trees. Its antipode is given by the recursive formula 

S*(t) = -t- (Id ★£*)(*) = -t - (S**Id)(£), (42) 

for £ 7^ i , and *S* ( i ) = i . 

To define the convolution of x(t) and y(i), we needed 
the condition x(\) = y(\) = 0. When this condition is not 
satisfied, we have two solutions. The first solution is to 
isolate the root, so that 

X(\)Y(\)=x( l )y( l ) + (X(\)-x( l )Mi) 

+x( l )(Y(X) - y(,) + (X(X) - xW)(Y(\) - y(,)) 

= x( l )y( l )+x( l ) * W V® 
\t\>o 

+ ^Al t lx(%(,)+^Al*l(a ; *y)(t). 

|t|>0 |t|>0 

The second solution is to use the coproduct A p . Thus, 
we define the convolution V by: 

(xVy)(t) = ^] x(uj)y(vj), where A p (t) = ^^Uj © Vj. 

3 3 

With this alternative convolution, the equality 

X{X)Y{X) = Y,^K^y){t) (43) 

t 

is satisfied even if x or y is not zero on the root. 

In our final formulas, we prefer to use the convolution * 
because its ensures the recursivity of the expressions (the 
trees in (x*y)(t) are strictly smaller than t). 

A last point of notation. If cp and ip depend on other 
arguments, we leave them inside (p and ijj. For example 

n(t) 

(ip(q)*i/>(q))(t) = ^2^{uj-q)^{vj-q). 

3 = 1 
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4.6 The self-energy 

As a first application of the convolution defined in the 
previous section, we introduce the tree-expansion for the 
electron self-energy. 

The bare electron self-energy S(q) is defined by: 

S-\q) = 7 • q ~ m - Z(q) = - V ef 'v»°(t; q), 



t 



where 



tl>°(r,q) = -('Y'q-m) and S(q) = ^ ef q), 

\t\>o 

so that ; q)ip°(\ ; q) = — 1. The pruning operator is 
used to define the expansion of the bare self-energy over 
trees in terms of the expansion of the bare electron Green 
function over trees: 

iP°(t)=i> o ( l )tp o (t)r(0+Tp O (i)(<P O *Tp°)(t), (44) 
<p\t) = ^(OVWO) + (¥>° *tf )(*)¥>°(')- (45) 

In terms of the antipode, ( p4|) can be rewritten 

V°(i) = -V ( 1 ;g)^ (g)o5,)(i)^ (i;g). (46) 

Similarly, for the renormalized self-energy, we have 

^°(t) = ; <?)(/(<?) o 5*) m°( , ; g). (47) 

We must give some detail concerning the meaning of ex- 
pressions like (cp°(q)oS*)(t). Because of its definition ((42|), 
the antipode S* acting on t generates a sum of products of 
trees. The action of (p°(q) on this sum is prolonged from 
its action on Y to an algebra homomorphism over T(Y). 
In other words 

/ (<?)(*! +£ 2 ) =<p°(t 1 iq) + (p°(t2iq), 
<p°(q)(\t)=\<p°(t;q). 

For a product of trees, we do not want to simply multiply 
two Feynman diagrams for the electron propagator, we 
must cancel one of the free propators between them. Thus, 
the product is 

f°(q)(tit2) = V(«i;9)^°(';«V (*2;«). (48) 

This operation becomes clear if one tries it on some ex- 
amples given in appendix 3. Since the (matrix) product 
on the right-hand side of Eq.(^) is not commutative, the 
algebra product on trees is not commutative either. 

In the presence of an external field A, the definition 
( |l5| ) of S°(z, y\ A) gives, after inversion and Fourier trans- 
form, 

il>°(\;q;A) = -(j-q-m-ej.A(q)) (49) 

Thus we obtain, at A = 

rl>°(r,q) = -(<y a q a -m), 
^(i;g;A,p) = 7 A , 
^ n (i;g;{A,p}i, n ) =0 for n > 1. 



The components of ijP(t;q) for the other trees t are 
obtained by using the chain rule for the functional deriva- 
tive of (Q) with respect to eoA^.fjpi), taken at A = 0. 
For n = 1, this gives the same res ult as in sect. 6.4 of [ [To| . 

Finally it is shown in section 13. 3| that the bare self- 
energy can be calculated from the recursive equation 



ip°(t;q)=i 



d 4 P 
(2tt)4 



7 Vaa' (ti ; P)9(t r ; q-p;X',p), (50) 



where 



g(t r ; q - p; A» = -(/(<? - p; \' ,p)V^°(q))(t r ) 
= &°(q)V^(q-p;\',p))(t r ). 

4.7 The higher components (f n (t) 

For a complete recursive solution of the renormalized prop- 
agators, we must define the higher components flHi) and 
(p^ u (t). As for the bare propagators (see Ref. (Uj), they 
are defined as the functional derivative with respect to an 
external electromagnetic field. 

Let us be more accurate concerning this external field. 
As noticed by Bogoliubov and Shirkov, Eq.([l6|) is not 
the Schwinger equation for QED with an external electro- 
magnetic field, since the latter involves tadpole diagrams 
which are absent from Eq.([L6|). However, Eq.(|l6|) is the 
Schwinger equation for a renormalizable theory (i.e. QED 
without tadpoles), and Dyson's relations (|8|-p!2|) still hold. 

In the real space, the bare and renormalized electron 
Green functions are expanded as 



S(x,y;A) 
S(x,y;A) 



t 



In these expressions, we do not distinguish between A and 
A because A comes always multiplied by eo and eoA° = 
eA. On the root, we have ip° ( i ; x, y\ A) = S°(x,y; A). 
The higher components of <p°(t; x, y\ A) must satisfy 

(p n (t;x,y;{\,z} 1:n ;A) = 



eSAx(z) 



(t;x,y;\,z,{\,z} lyTl ;A). 



whe re the notation (f 7 ^ 1 and {A, z}i^ n is defined in section 

Of 

Since our purpose is QED without external field, A is 
just used to take functional derivatives, and the higher 
components we actually need are 

(p n (t; x, y; {A, z} lyU ) = (p n (t; x, y; {A, z} hn ; A) 

for A = 0. 

At A = 0, the theory becomes translational invariant, 
and a Fourier transform gives us 

(p n (t;q;{\,p} hn ) = ip™ +1 (t;q; \,p,{\,p} hn ). 



eSA x (p) 
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In the recursive equations for (p(t) we meet products 
of propagators such as (p°(ti;q)(p°(t2;q). In the real space, 
this gives a convolution of cp°(ti; x, y) and (p°(t2'<, y). We 
take the functional derivative of this convolution with re- 
spect to A(z), we Fourier transform the result and we 
obtain, at A = 

0°(ti; q)<p°(t 2 l q)) = ^faq)? 1 (t 2 ;q; \,p) 



e5A x (p) 



+ip 1 (t 1 ;q; \,p)tp° (t 2 ; q + p). 



This expression satisfies energy-momentum conservation. 

To take the functional derivatives of the recursive equa- 
tions for renormalized quantities, we need the indepen- 
dence of the renormalization constants with respect to the 
external field. There are various ways to prove this. For in- 
stance, the differential form of the Ward identity (Eq.(21) 
in Ref.Jnj is 



d(p n (t',q]{\,p} hn ) 
dqv 



— _^ n + 1 



(/4 +i (£;g;/i,0,{A,p}i, n ). 



The Ward identity is also valid for the renormalized elec- 
tron propagator |57| , so that 



d(p n (t',q',{\,p} ljn ) 
dqv 



(/4 +1 (^;/i,0,{A,p}i, n ). 



From the definitions ( [46|) and (p7| ) of the bare and renor- 
malized self-energies we obtain 



dip n (t; q; {X,p}i, n ) 

d<p n (t;q;{\,p} hn ) 
dqv 



iPl +L (t;q;fi,0,{\,p} 1)n ), 



-^ +1 (t;q;p,0,{X,p} hn ). 



Now we start from the relation between the renormal- 
ized and bare self-energies, for instance 

V;°(Y;<7) = V°(Y;<?) - C 2 (Y)( 7 - q - m) - < m (Y). (51) 

On the one hand, we take the derivative of Eq.(|5l|) with 
respect to q^ and use the Ward identities (and the fact 
that C2OO and £ m (Y) do not depend on q) to obtain 

- ^(Y; q- 0) = -^(Y; q; /i, 0) - C 2 (Y) 7 M . (52) 

On the other hand, we take the functional derivative 
of Eq.(|5l|) at A = and we obtain 

^\y;q;H,p)=i; 1 (r,q;p,p) + C2 (Y) 7 M - C2 00 (7 • « - (8?) 

-Cm 00, (54) 

where C2OO an d Cm 00 denote the derivative of C2OO and 
Cm(T) with respect to A^(p) at A = 0. The term C 2 (Y)7^ 
comes from the functional derivative of Eq. ( p9| ) . 

If we take the value p = in Eq.(p3) and compare 
with Eq.(||) we obtain C2 00 = and Q^f) = 0. Further 



differentiation shows that Q (Y) = and C™ 00 = 0. 
We can apply this proof to any tree t, once the subdiver- 
gences have been subtracted. A similar proof can be given 



for (3(f), using Fury's theorem instead of Ward identities. 
This proof assumes that the renormalization conditions 
do not depend on A (e.g. minimal subtraction). 

More physically, S(x,y; A) has the same singular struc- 
ture at x = y as S°(x,y), except for logarithmic terms 
that are integrable. Thus, the renormalization constants 
are determined by S°(x,y). 

This method of functional derivatives avoids the usual 
renormalization of vertex diagrams. Renormalizing prop- 
agators is sufficient. 



5 The renormalized electron propagator 

In this section, we show that the recursive equation for 
the electron propagator is 



<p°(t;q) = p(t)<p°(r,q)-( m (t)<p°(r,q) 2 
V(';«)(Cm *£"(«))(«) 



+V(';<?) 



d 4 p 

W) 4 



7 A (pl x , {ti;p)(p 1 {t r ;q—p; A',p),(55) 



where p(t) = — C2(*) - (p*C2)(*) = C2 S*(t), starting at 
p(Y) = -C 2 (Y). 

It is natural to define a new quantity a(t] q) by a(t; q) = 
for t = 1 and, for t = ti\/t r , 

■j—y 1 X ¥\\> (U ;p)(p 1 (t r ;q-p;\',p). 

Then we consider Eq.(p~7|) and in the integral over p, we 
expand the photon propagator over trees t\ using Eq.(|24|) 
and the electron propagator over trees t r using Eq.(p3[). 
We recognize a sum of a(U V t r ) and the integral becomes 



ie'S^q) 



d A p 



J X Dxx>(p) 



SS(g-p) 
eSA x >(p) 



^2<x(t]q). 



In the other terms of Eq.([l7|), we expand S(q) and 
renormalization constants over trees using Eqs.(^), (|3l|), 
"~2|) and (33). We replace products by convolutions V ac- 



cording to Eq. (|43|) in all expressions except the integral 
and we obtain. 

5>i*V( g )vc2)w = ¥>°(i;«) + E elt| («(9) v C2)W 

t t 

-£eV(';«)(CmV$>°(g)VC2)(t). 

t 

The bold step is now to identify the terms corresponding 
to a given tree t. This yields 

($>°(«)VC2)(t) = ¥>°(i; 9 )6(t) + (a(q)V( 2 )(t) 

-f°(r,q)(C m V^(q)VC2)(t), (56) 

where e(t) = 1 if t = i and e(t) = otherwise. 

To simplify this expression, we follow Kreimer []58| and 
compute (^°(^)VC2VC2 oS*)(t). The basic property of the 
antipode is IdX7S* — ie, therefore 



C2VC2 o S* = C 2 (/dVS*) = C 2 e = e. 



(57) 
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The associativity of V is here crucial. On the one hand, 
(£°te)V(C 2 VC2 o £*))(*) = <p°(q)(t) according to @. On 
the other hand, we calculate ((<j? (<z)VC2)VC2°5'*)) , where 
we replace <g°(g)VC2 by the right-hand side of Eq.( f56|) . 
Equation ( |57| ) gives us 

((^ (9)VC 2 )VC2 °SJ) = C2 o5,(%°(,i9) 

V(';9)(C m V/( ( /))(t)+a(t). 

From the associativity of V and the definition of a(t; q) 
we obtain our final recursive equation ( |55| ) for the electron 
propagator. 

The recursive equation is completed by the equation 
for the higher components of (p°(t; q). If we take the func- 
tional derivative of Eq.(p5|), and make the same simplifi- 
cation as in Ref . Q ) , we obtain 

^(t;q;{\,v}i,n)=^{r,q) 1 x ^ n -\t;q ] {\,p} 2 ^) 

-p(t)<p°(\;q)6 nfi - Cm(t)(p°(i;q)<fi n (i;q;{X,p}i, n ) 
-<P°(r,q){C m *^{q;{\,p}i, n ))(t) 

W(';?)2 / 77^47 A ^AA'(^;p;{A,p}i,fe) 

fe=0 ^ ' 

x^l~ k+1 (t r ; q-p; X',p, {X,p} k+hn ). (58) 

Another useful formula can be obtained by defining 

f°(t;q) = &°(q)V( 2 )(t) (59) 
= <p° (t ; q) + &° (q) * ( 2 ) (t) + cp° ( , ; q) ( 2 (t) , 
P(t;q;\',p) = {(p\q- A',p)VC 2 )(t). (60) 

With this notation Eq.(j56|) is rewritten 

?) = -Cm i ; <z) 2 - ¥>°0 ; «) * /°(«))(*) 



+iip°(r,q) 



d 4 p 



rvlx(tr,p)f 1 (t r ;q-p;\',p). (61) 



i V t r . Now we show that the recursive equation for the 
renormalized photon term ^ A (i V t r \q) is 

Vt;«) = -Cs(i Vt)^(i; 9 ) 
d 4 p 



-^VaO;?) 



tr^^^A',-?)]^,;^) 



where J 1 was defined in the previous section. 

To prove this, we start from several remarks: we have 
Vi r ) = is ^ i and 

^('Vt r ;9) = -< A (i;g)£ 0AV (i Vr^Lt ';<?)• 

Because of this close analogy between photon propa- 
gator and vacuum polarization, we shall rewrite (^l|) as 

5> 2 l^(, vt) = -^.;g)77 AV (^0;9) 



-(Z 3 -l)^(i;«) 



^2«e <^ A (i;g) 



d 4 p 



tr[7 



<J5(p) 



(2tt) 4 l ' eSAy(—q) 



|vvi/(';«)- 



We rewrite this expression to isolate the root compo- 
nents: 



Vt) = -(Z3-1)^(,; ? ) 



(2tt)4 l' eM v (-g) J 

x^v(';?)- 



We expand all quantities over trees, using 
6S(p) 



eSAy(-q) 



J2e^\t;p;X',-q), 



The higher components f n (t;q) are obtained by func- 
tional derivative of (|6l|), as explained in sectio n |4.7| . The 
recursive equation for f n (t; q) is the same as Eq.(|58|), where 
if is replaced by / and the term p(t)cp°( i ; q)5 n ,o is sup- 
pressed. 



6 The renormalized photon propagator 

Bogoliubov and Shirkov j2lj have shown that the renor- 
malization of two Feynman diagrams linked by a single 
photon (or electron) line is obtained by an independent 
renormalization of each of the two subgraphs. In our lan- 
guage, this means that the renormalized form of ( ^9|) is 

<P%,(ti Vt P ;«) = £° A (' V ■t r ;q)il> oxx ' '(i;«)^(*«;9)(62) 

Therefore, all trees for the photon propagator can be 
renormalized once we have renormalized the special trees 



and we multiply through the pruning operator. Then we 
identify the terms corresponding to the same tree and we 
obtain 

Vt;g) = -CsO Vt)^(i; g )+C2W^(Y;g) 

— -^tr [ 7 V * C2)(«;p; A', -<?)] ; g). 

From the definition ( |60|) for J 1 , we can rewrite this 
expression as our recursive equation (|63|). 

The higher components are obtained very simply by 
taking the functional derivative of Eq.(|6^). Since (pj^ u ( i ; <?) 
is independent of the external field, we obtain 

^0 V£;g;{A,p}i, n ) = -i^ A (i;g) 
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For the other trees, we use 

n 

<f%v(U vfr;g;{A,p}i,„) = J^*a(' vi r ;g;{A,p}i, fe ) 

fc=0 



xV 0AV (i;g)^- fe (t I ;«;{A,p} fe+ i,„). 



6.1 Properties of renormalized photon propagator 

From (|6l| ) and ( |63|) we can deduce that the renormalized 
photon propagator does not depend on any C2^)- In fact, 
we shall prove that f°(t;q), J 1 (£;<?; A,p) and (p^ u (t;q) do 
not depend on any (2(f)- To do this, we reintroduce a 
non-zero external field A. The property is clearly true for 
t = 1. If it is true for all trees with \t\ < N, let us take 
a tree with \t\ = N. Because of (|6l|), f°(t;q) does not 
depend on any (2^)- Since / 1 (t;^;A,p) is obtained by a 
functional derivative of g) with respect to eA, it does 
not depend on any either (eA does not depend on 

any C2 (£'))• If t = 1 Vt r , because of (|63|), Vt r ;q) does 

not depend on any £2^) since none of the terms on the 
right hand side do. Finally, if t is not of the form t = 1 V£ r , 
it is of the form t = UV t r , and (p^ifi V t r \q) is obtained 
from <£>^(i Vt r ;q) and (p^ v (ti] g), which do not depend on 
anyC 2 (i')- 

With the same reasoning, we see that f°(t;q) and 
(f°^{t;q) are independent of the gauge parameter £ for 
f ^ 1. 



7 Electron self-energy 

To calculate the electron self-energy, we start from ( ^5|) 
that we rewrite 

<p° (t;q)= p(t)tp° (\;q) — Cm(t)f°{r,q) 2 + a(t;q) 

V(>;«)(C ro */(«))(*)■ 

The self-e nerg y is obtained by introducing the last 
equation into ((44|). This gives us 

^°(t; 9) = -pm°( 1 ; <?) - Cm(t) + V>°( 1 ; «)<*(*; <?)V>°( ' ; <?) 

-( P *4,°( q ))(t) + (t m *^(q))(t)i> (r,q) 

+<P%;q)Km* ^° («))(*) + (Cm * <P°(q) * ^° («))(*) 
+r(r,q)(a(q)*4,°(q))(t). 

If we factorize Cm and use (|45|), the expression reduces to 

$°(t; g) = - P m°( , ; <?) - C™ W + V>°( ' ; «)<*(*; <?)V>°( ' 5 <?) 
-(p* ^°(«))(t) + ; «)(<*(«) * 4>°(q)) (t). 

From the definition of a(t; q) and the result of section 



13. 3| , we obtain 

^ (*; q) = -pm° ( i ; q) - Cm (*) - (p * ^° (9) ) (*) 

/d 4 p 



where 

g(i r ; g - p; A',p) = -(^(g - p; A',p)V^°(<?))(t r ) (64) 
= (<fP{q)Vp{q-p;\',p)){t r ). (65) 

It can be shown recursively that the only term pro- 
portional to ip°(\',q) in —p(t)ip (\;q) — (p * ip (q))(t) is 
(2m°(r,q). 



8 Renormalization conditions 

Renormalization conditions are the conditions used to give 
a unique value to the renormalized quantities. As noticed 
by Itzykson and Zuber there is some freedom in the 
choice of these conditions. However, it is necessary that 
the conditions are satisfied for the root. In other words, 
the renormalization conditions can be non-zero only for 
the root. The physical meaning of these renormalization 
conditions can be found in any textbook on quantum field 
theory. In particular, the relation between the renormal- 
ization conditions and the introduction of experimental 
uantities into the theory is discussed at length in Refs. 

For instance, mass shell renormalization conditions for 
QED Q p. 413 could be translated into (for t ^ i ) 

$ (t',q)\<y» qti =m = 0, 
^ (£;<?; A, 0)| 7 /^= m = 0, 
^0;<?)|g=o = 0. 

The strange prescription 7 /x # /x = m means that the 
quantities $°(t] q) and q\ A, 0) must be multiplied by 
u(q) on the right where u(q) is a solution of the Dirac 
equation (j^q^ — m)u(q) = 0. This amounts to replacing 
q 2 by m 2 and 7 M g M by m in the analytic expressions for 
$°(t; q) and q; A, 0), when they are available. 

For the photon, we define Co(q 2 ) by 

[ D ~\M = (9Ag A1 -9 2 5Ap)(l + ^(9 2 ))-^.(66) 
This Co(q 2 ) is the same as ujR^q 2 ) in Itzykson and Zuber 

i- 

For the photon, we know that i/j® u (t;q) is transverse. 
Therefore, we can define uj{\ V t;q) by 

V^O V^)=^(,;^(i Vt-q). 

Notice that Co(\ V t\q) = Cu(q), where uj(q) was defined 
in (|66|). The renormalization condition on Cd(\ Wt;q) re- 
places a more complicated renormalization condition on 
ip^vfaq) (involving the third derivative of ijj {t\q) with 
respect to q). 

To show how this works in practice, we use the identity 
Vt;«) = -w(i Vi;g 2 )^(,;g) 
Now we assume that we have calculated the integral 



Hq) = ~g XX ' 



d 4 p 

W) 4 



tr[ 7 x f\t;p;X',-q)]. 
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This integral is logarithmically divergent, and we must 
determine how to remove its divergence. We start from 
Eq.(|63|), that we multiply by g^ v and sum over \i and v. 
Using the fact that the integral over p is transverse, we 
obtain 

3w(i Vt;q 2 ) =3Ca(i Vt) + %). 

Since the renormalization condition is 0(\ V t\ 0) = 0, we 
just take (3(1 V t) = —h(0)/3. This satisfies the renor- 
malization condition and this determines ("3(1 Vt). For the 
electron propagator, the renormalization conditions deter- 
mine ( 2 (t) and (m(t). 

In practice, other renormalization conditions are used, 
such as minimal subtraction. Once quantities have been 
made convergent by a renormalization condition, it is al- 
ways possible to translate the results into different renor- 
malization conditions, using the Hopf structure of renor- 
malization ]58[ . The fact that the renormalization condi- 
tions can be composed and that the result does not depend 
on which condition is used first is a direct consequence of 
the coassociativity of the Hopf algebra of renormalization 



9 Massive quenched QED 

Quenched QED is QED without vacuum insertion graphs. 
Quenched QED has been recently advocated and discussed 
by Broadhurst and coll. (59||6(J. 

The Feynman diagrams describing the electron propa- 
gator of quenched QED have no fermion loop. Within the 
present approach, all Feynman diagrams of a given order 
for quenched QED are summed into one tree (/?°(Y n ). The 
trees Y n are defined recursively by Y° = 1, Y n = 1 V Y n_1 . 
In words, they are the trees without left branches. The 
Feynman diagrams corresponding to (/?°(Y n ) for n < 3 are 
given in appendix 3. 

The bare and renormalized electron propagators for 
quenched QED are written 

00 00 
S«(g) = 5>fa>°(Y";«) and S%) = £ eV(Y"; <?)• 



n=0 



n=0 



The Schwinger equation for the renormalized electron 
propagator of quenched QED is 

S Q (q) = S°(q) - (Z 2 - l)SQ(q) - 5mZ 2 S° (<?) 

■He Z 2 S (q) J ^ 7 D xx ,(p) . (67) 

The relation we want to show is 



n—l n—a 



IV j_ i v i x -| 

gP^q) = <p\r-q) + E E k^ n ~ a ) 



a=0 k=0 



v 1 u n 9y(g) VYB) 



(68) 



where 

a = p = -C2 - P*C2, 

Oil = ~Cm ~ Cm *«0, 

afe = -Cm *afc-i if fe > 2. 

Notice that o^(Y n ) = if k > n, because Y n cannot be 
split into more than n part by the pruning operator P. 

To prove this, we start from the Schwinger equation 
( |f37|) and transform it into a recurrence relation: 

<p°(t;q) =p(t)v°(r,q)-( m (t)<p°(r,q) 2 

-p°(';«)(Cm */(«))(*) 
d 4 p 



-iS°(q) 



(2n) 



n X D° xxr (p)^(t r ;q- P ;X',p), (69) 



where p(t) = -C2W - (p*C2)(*) = C2 starting at 

P(?) = -C2(Y). 



Using the definitions given in section 4.7, we can show 



that 



d 



edA\(p) dm 
and that 



<p°(r;q) 



d 



dm edA\(p) 



V\^\q), 



^(^q)=^yq)^^pq) 

+i(f0{r ' q) J ^^^^^q-p-A'iP). 

From these identities, Eq.(^) follows easily by recursion. 

The Feynman diagrams describing the photon propa- 
gator of quenched QED have a single electron loop, as can 
be seen in the diagrams for ^ M (Y n ) in appendix 3. The 
bare and renormalized photon propagators are written 

oo 
n=0 

CO 

^(9) = E e VA,(Y n ; < /). 



n=0 

The relation we want to show is 

n—l n—k 



Vl^ n ;q) = Vl^ n ;q) + E E h.^ n ~ a) 



d^Sf^q} 



k=0 a=l 



(70) 



-C3(Y")^(,;g), 
where the scalars /3k are defined recursively 

Pi = "Cm, 

Pk = -Cm*Pk-i k>2, 



with /fe(Y n ) = if k > n. This can be written more com- 
pactly: 



fl^g) =exp[-C„ 



_d_ 

dm 



WU^q)-Cs(r n )^(r,q). 
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To prove Eq.@, we first use our previous result Eq.(p8 
to show that 



n n—k 



f°(r;q) = ¥>°(Y";9) + £E ^0^"°) 



1 „ ( ^_ a ,d k ^{^-q) 



k=l a=0 



Then we introduce this expression into Eq.(p3j), and 
the result follows from the relation 



Qk 



I (^^^(Y-^A',-,)]^,^,;,). (71) 



10 Hopf algebra for massless QED 

In this section, we determine a coproduct from the re- 
cursive equations f55|), ( |6l| ) and ( |63| ) and the product law 
( |62| ) for massless QED. The case of massless QED is much 
simpler because the mass is not renormalized. Thus, for 
all trees t, ( m (t) = 0. 

This coproduct determines the renormalized propaga- 
tors as a function of the unrenormalized ones. In this sec- 
tion, it will be useful to dinstinguish the electron and pho- 
ton trees by the color of the root. A tree with a black root 
is written t* and represents an electron propagator, a tree 
with a white root is written t° and represents a photon 
propagator. In a tree t\ V t r , t\ is white and t r is black. 
There are now two graftings operators V and V , so that 
tiVt r is a black tree and ti\/t r a white one. 

Using a variation of Sweedler's notation, we write 



4t° = E*(l)® f (2). 
At° 

At* = ^^(1)^(1) ^(2)5 



Af 



F (t')= E *(1)®*(2)- 
F(t') 



(72) 
(73) 
(74) 



These equations mean that the coproduct of t° generates a 
sum of tensor products with one white tree on the left and 
one white tree on the right, the coproduct of P generates a 
sum of tensor products with one black tree and one white 
tree on the left and one black tree on the right, finally the 
coproduct F(t) generates a sum of tensor products with 
one white tree on the left and one black tree on the right. 
These trees can eventually be the root, which is the unit 
element of the algebra (the root is neither white nor black, 
or both, as you wish). 

To avoid products of white trees in Eqs.(f72|), ( |73| ) and 
Qj_we took advantage of the fact that, according to 
Eq. (p29|) , the cp^ of a white tree ti\/t r can be written as 
a product of (p^u (iV£ r ) by (p^ u (ti). From Eq.(p2l), we also 
know that this property is compatible with renormaliza- 
tion. Therefore, we translate this property into an inner 
product over white trees. The product of two white trees 
s°/t° (read "s over t") is defined recursively by s/\ — 



s and s/(ti\/t r ) = (s/ti)Vt r . In particular ti/(\Vt r ) = 
tiVt r , which is what we need. Surprisingly, this product 
has been used independentl y by Loday and Ronco in a 
completely different context |6l| . 

The coproduct A acting on white and black trees is 
defined by the recursive equations 



A(\Vt) = (iVt) <8> i + XXi) ( lVt (2))' 



(75) 



F(t) 



A(tiVt r ) = (t l Vt r )®l+ (*?(l)Ar(l))£(l) 

At u At r 

^(^(2)V^ (2) ), (76) 
F(t l Vt r )= (%)Ar(i))^(^°(2)V^ (2) ), (77) 

Ati,F(t r ) 

with the initial values A\ = i i and F(\) = i i, and 
with the compatibility of the "over" product with renor- 
malization: A(s°/t°) = At°As°. In particular, 

AtiVU = A(\Vt r )Ati. 

With this notation, we can now write the coproduct of 
a general white tree 



1(2) 



+ E (^(l)Ar(l))^(^(2)V^ 2) ). (78) 

At U F{t r ) 

These preliminaries enable us to write the relation be- 
tween renormalized and unrenormalized propagators as 



^U*°;9) = Ec(*(i))^(*(2);<z), (79) 

At° 

<fP(f;q) =YC(t° (1) )C(t' 1) )<p (t' (2 y,q), (80) 

Af 

mq) = J2((t° (1) )f°(f (2) ). (81) 



Af 



The general counterterm ( is a scalar over black and 
white trees defined by ("(0 = 1 and 



c(s°/n = c(s°)c(t o ), 

C(.Vt) = — Ca(i V*). 



(82) 



Eqs. 



In particular ((tiVt r ) — — ^3 ( 1 v U)((ti)- We recall that 
p(t) = C2 ° S*(i). Equations and (80) are given in 
expanded form in appendix 2 for trees up to order 3. 
We prove this recursively. From the list of appendix 2, 
and ( |8Q| ) are satisfied for all trees up to order 3. 
The same can be checked for Eq.(|8T|). Assume that they 
are satisfied up for trees with 2N — 1 vertices. Take a tree 
with 2A/~ + 1 vertices. Take first iVt, then use Eq.© for 
f 1 in Eq.(||). This yields 

Vt;g) = -C 3 (' Vt)(pl„(r,q) 

+ Z)C(*(i))^(i vt (2) ). 

F{t) 
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This is Eq.@ for the coproduct defined by Eq.@. If 
we take now ti\/t r , where t\ and t r have less than 2N + 
1 vertices, we can expand (p^ v {\ V t r ) and (p^iti) over 
unrenormalized terms. Then, using Eq.(|62|) we find 

<pl„(tiVtr,q) = (-(C 3 (i Vt r )^(,;g) 

+ C(*z(i)Ar(i))^(^(2) vt r(2 );g). 

F{t r )AU 

This is Eq.(f79[) with the coproduct defined in Eq.([78|) and 

the correct Q. 

For the coproduct acting on electron trees, we start 
from the recursive equation ( |55| ) and we use the expansion 
© for (flvfcq) and © for <p°{t r ;q). This gives us 

<p°(t r Vti',q) =p(tiVt r )(p°(r,q) 

+ E C(*i(l)Ar(l))C(*; ( l))^°(*i(2) Vtr(2)i?). 
F(t r )Att 

The first term p(t / Vt r ) is consistent with V t r ) as 
defined in Eq.(p2|). And the other terms are consistent 
with Eq.(pO|) using the coproduct Eq.(|76|). 

Exactly the same substitution leads to Eq. (|8l|) using 
the coproduct Eq.(^). 

In Ref . ]6^| , it will be shown that A is coassociative 
and defines a Hopf algebra over the two-coloured planar 
binary trees. 



at each order of the perturbation theory. Each order of 
perturbation theory of quenched QED is renormalized by 
the standard renormalization of Feynman diagrams, so we 
know that our (^°(Y n ; q) and ^-^(t; q) are finite, since they 
are uniquely fixed by the Schwinger equations. The other 
trees are obtained from some Y n by insertion of renormal- 
ized photon propagators. Since these insertions do not in- 
troduce new divergences, and particularly no new overlap- 
ping divergence, we can conclude that all trees are finite. 
The point that is still not solved is the insertion of higher- 
order photon Green functions, such as the photon-photon 
scattering diagrams. Since these diagrams are finite, they 
probably do not spoil the convergence, but we could not 
prove this generally. 
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13 Appendix 1 

This appendix contains some proofs. 



11 Conclusion 

The method of Schwinger equations has a number of ad- 
vantages: operator- valued distributions are avoided, as well 
as indefinite norms and many of the difficult mathematical 
concepts of quantum field theory. As compared to Feyn- 
man diagrams, the method of planar binary trees is more 
compact and does not require symmetry factors. More- 
over, our treatment of renormalization does not require a 
special treatment of the so-called overlapping divergences. 
Such a special treatment is even necessary for Kreimer's 
original method of renormalization by rooted trees Jp8| . 
The present work was much inspired by Refs.|j||f7[ 
but it was carried out independently of the recent and 
fascinating results by Connes and Kreimer Still, our 

Hopf algebra can probably be modified to fit into their 
general framework. We hope to explore this connection in 
a forthcoming publication. 

We conclude this paper with a word of caution. We 
have not actually proved that each tree is renormalized by 
the Hopf algebra. In other words, we have not proved that 
all <p°(t; q) and (fx^t; q) are finite. However, we can give an 
argument that can eventually lead to a proof of finiteness. 
In the case of quenched QED, there is a single tree (Y n ) 



13.1 Proof of @j (the pruning operator is 
coassociative) 

If t = i we have P(\) = 0, so the identity ( |4l| ) holds. So, 
suppose that t ^ i . The reason why the identity ( f4l| ) holds 
for t is that applying the recursive definition of P, on the 
successive right grafters of t, at each step both sides of (|4l| ) 
coincide on the terms which do not involve P(t') for the 
last right grafter t' considered. Of course, when we finally 
meet a right grafter t' such that P(t') = we obtain the 
equality (jll|). We develop this idea formally. 

Any tree t ^ i can be written in a unique way as 

t = h v(t 2 V(... V(t n V i)...)), 

for some n < \t\ + 1. In fact, it suffices to decompose the 
tree t into its left and right grafting trees, then to decom- 
pose successively the right trees as graftings of two new 
trees and pick up all their left sides, t\ := fj, t 2 := {t r )h 
h '•= ((t r )r)i and so on, until we meet an undecomposable 
right side (...((£ r ) r )...) r = i. Since \t\ = \t ± \ + \t 2 \ + ... + 
\t n \ + n — 1 and |^| > for all i = 1, n, the procedure 
must finish for an n < \t\ + 1. 
Since P(t n V i) = 0, we have 
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P(t n -1 V(£ n V i)) = t n -l V I ®t n V I, 
P(*n-2 V (t n _i V (t n V I ))) = t n - 2 V I t n -l V (t n V I ) 

+ t n -2 V (t n _i V l) ®t n V I 

and 

P(t n - 3 V (t n _ 2 V (t n _i V (t n V i)))) = 

t n -3 V I ®t n -2 V (t n -l V (t n V i)) 
+*n-3 V (t n _ 2 V I ) (8) £ n -l V (t n V ! ) 
+*n-3 V (t n -2 V (t n _i V I ))®t n V I 

Thus, for the tree t = t\ V (t<i V (... V (t n V i )...))> we obtain 

n-l 

= v(t 2 v(...v(tiV i)...)) 



iV(?) = E e oV^) 



e(^) + E e o |tll+2|£2l+2 ^(^Vt 2 ;,) 



tit 2 



P. 



^(';^) + E e o 



2|ti|+2|t 2 |+2 o 



tlt 2 



2=1 



H+l 



V (t <+ 2 V (... V (t n V l)...))- 



Hence 



(P (g) irf) o P(t) = 

n-l 

= ^P(i 1 V(i 2 V(...V(i,V ,)...))) 

j=l 

V (ti+2 V (... V (tn V I )...)) 
n-l J-l 

= £X>V(...(tiVi)...) 

J=l »=1 

®t i+ i V (...(tj V i)...) ®fj+i v (...(*„ V i)...) 
53 tiV(...(tiVi)...) 

l<i<j'<n— 1 

(8)t i+ i V (...(*, V i)...) V (...(*„ V [)...), 
and similarly 
(id^P)oP(t) = 

n-l 

= E^i v fe v(...v(t fc V i)...)) 

(g)P(t fc +i V(* fe+ 2 V(... V(t„ V i )...))) 

n— 1 n— 1 

= EE *iv(...(t fc vi)...) 

fc=i i=k+i 

®tfc+i V (...(t, V i)...)®ti+iV(...(i„V i)...) 

51 *i(...v(t fc v i)...) 

l<fe<J<n-l 

®t fc +l V (...(*, V i)...)®i i+ iV(...(i„V i)...). 



Hence the identity (41) holds for any tree. 



13.2 Proof of @ and (§) 

We use © to prove M). 



xV 0AA '(i;«)^ v (ti;«) 
= ^( l ;9) + E e o t2l+2 ^A( l vt 2 ; (/ ) 

x^ 0AA '(,;g)I> Vl ,(g). 

We multiply the last equation by D~ 1 (q) on the right 
and by Z}° on the left. This gives us 

[D°-\M = [D-\M) + T.^ +2 [^-\ x { q ) 

t 2 

^ XX '(^t 2 ;q)^(r,q). 

Since (p oxx ( i Vt2; q) is transverse, we can replace [D° 1 } /J/ \(q) 
by ^ A (i;g) in the above equation. 



Using |28| and the definition of U^ v (q) given in ( |25|) we 
obtain 

*W«) = "E e o" 2l+2 ^x(r,q)v° XX '(< Vt 2 ;q)4> x ,„(r, q ). 

t 2 

We can rewrite this last equation as 

n,M = E ef = ^ ef l+2 <„0 Vt;g), 

|t|>o * 

with V t 2 ; g) = iffi ^ i and 

i V t; g ) = -^ A ( 1 5 <zV° AA ' ( i V t; g)^,„( , ; g). (83) 

Finally, we use (|37|) with n = and £1 = i 



Vt 2 ;q) = -iDl x {q) f ^tr[ 7 ^\t 2 ;p;X',-q)] 

We multiply this equation by D° 1 (q) on the left and 
on the right, and we use the fact that V t2',q) is 

transverse to get, from (|83|), the relation 

^"(,Vt;g)=i J ^tr[ 7 V(^, "<?)]• 
If we sum over trees t, the last equation becomes 



d4 P ...r a ss (?) 



-tr[ 7 
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13.3 Proof of © 

From the definition of self-energy ?/? (Eq.(|44|)), we use the 
definition of the convolution (Eq.(p9|)) and of the pruning 
operator (Eq.(f40|)) to write, for a tree t = £j V t r 

= tp°{>)<p°M%) + V> (.)(/*V )W 

n{t r ) 

+ £>%y ) (tjVu i ty (t; i ), 

2=1 

where ui and V{ are the trees obtained by pruning t r (i.e. 

The last equation will be transformed by using the 
recursive definition of cp° for the trees t, ti V i and V^. 
For n = 0, Eq.(|35|) becomes 

<P°{t;q)=i J ■^S\q) 1 \l y {t l] p) V \t r - q -p-\' ) p). 
Therefore, we obtain 

^(tr'^q-p; X',p)i/j (\;q) + (/? 1 (i;<? -p; A',p)^°(£ r ; q) 

n(t r ) 

+ ^2 ^(uiiq-piX'iPj^iviiq) . 

2=1 

This can also be written 

/d 4 p 
^^tVaa' (U;p)g(t r ; q-m A',p), 

where 

g(t r ;q-p;X f ,p) = -^(t r ; g - p; A / ,p)'0°( i ; g) 

-^(ijg-p; \',p)ijp(t r ;q). 

It will be useful to transform g{t r ;q — p; X\p). To do 
that, we rewrite the equation for ^(t; #5 A,p) given at the 
end of paragraph 6.4 in Ref.]l0|, so that it gives 

<p°( i ; ^ (t; g; A,p) = <^°( i ; V(*; Q + ' ; 9 + p) 
-y>°(t; g) 7 A - ^ (t- q- A, ( i ; q + p) 

n(t) 

+ X! K ( 1 ' V ^ + ^)^° (v» ; g + p) 

2 = 1 

-yPiuiiq^tynqiXip) - ^ 1 (^;g;A,p)^ (^;g + p) . 

Now, we replace <^°(£; <7+p) by its value given from Eq.©, 
we add g(t;q; A,p) on both sides and we reorder a bit. 

<j(t; (z; A,p) = <p°( i ; g)^; g; A,p) + g) 7 A 

V(';^V(^ + ^°(';^p) 

-^(i;^; A,p)^°(i;g+p) 

n(t) 

+ ^ ^ ( Ui 5 ^ + P ^ ( Vi 5 ^ A ' * 

2 = 1 



Finally, we note that i ; g; A,p) = ; q)j X (p°(t; q+p) 
and we obtain 

g(t;q;\p) = ^O;^ 1 ^; A,p) + ^°(t;g)^ 1 (i;g; A,p) 

n(t) 



14 Appendix2 

In this appendix, we collect the relation between bare and 
renormalized photon and electron cp up to three loops. 

14.1 Electron Green function for massless QED 

14.1.1 One loop 

£(Y) = ^(Y)-C 2 (YMi). 

14.1.2 Two loops 

^ = ^a)-C3(YMY)-c 2 (^M'), 

<p(f) = p(f) - C2(Y)<^(Y) - C2(Y>(,) + C2(Y)V(0- 

14.1.3 Three loops 

= ¥>(>)- 2Cs(YM^Y) +Cs(Y)V(Y) - C^MO, 

^) = ^)-C3(fMY)-c 2 ^M0, 
^(Y) = ip(Y) - Cs(Y)^(r) - c 2 (YM > t) 

+C 2 (yR 3 (yMy) - C 2 (YM,) + C 2 ( s y)C 2 (y)^(0, 

^ = ^(T)-C 3 (YMr)-C 2 (>MY) 

-c 2 (YV(.)+c 2 ax 2 (YMo, 

= <p(f) - C 2 (YMY-) - C 2 (Y-MY) 
+C 2 (Y)V(Y) + [-(2(f) + 2C 2 (r)C 2 (Y) - C 2 (Y) 3 ]^(i). 

14.2 Photon Green function for massless QED 

14.2.1 One loop 

^M = ^M-(3(yW(')- 

14.2.2 Two loops 

^(^) = - 2C 3 (yVam(y) + CsMVa^O), 

^00 = ^(f)-C3(f)^(i)- 
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14.2.3 Three loops 



14.4.3 Three loops 



£am(^=<M>)-3C3(yW^) 

+3C 3 (T)Vap(Y)-C3(T)Vap(0, 
fxrtf) = <P\»Cf) - C 3 (y)^(y') - C 3 (t> A p(Y) 

+C3(f)C3(Y)^(0, 

<^(Y) = ^(Y) - C 3 (Y)^Ap(f ) - C 3 (f WW 

+C3(f)C 3 (Y)^(i), 
^(Y) = ^Am(Y) - C 3 (Y)^(f) - C 3 (Y)^(0, 
=^ v (f)- (3(^)^(0- 

14.3 Electron self-energy for massless QED 

14.3.1 One loop 



^(Y)=V(Y)+C 2 (Y)V(0- 



14.3.2 Two loops 

^(^)=v( s t)-C 3 (Y)V(Y) + C2( x Y)V'(0, 

^co = vco+c 2 (y)v(y)+c2W(i). 

14.3.3 Three loops 

<K>) = - 2C 3 (y)V(^y) + C 3 (y) 2 V(y) + C2(>)V'(0, 

^)=V(^)-C3(y-My) + C2^M0, 

Vi(Y) = V(Y) - C 3 (YMf) + C 2 ( s y)^(y) + C 2 (YM0, 

^(Y)=^(Y)-C 3 (YWr) + C 2 (Y)^( 1 ) 
+C 2 (Y)^( S Y)-C 2 (Y)C 3 (Y)V(Y), 
^(f) = ^(f) + C 2 (YMf) + C 2 (f)^(Y) + C 2 (^(i). 

14.4 Vacuum polarization for massless QED 

14.4.1 One loop 

^M=^(T)-C3(T)^(i). 

14.4.2 Two loops 

^m = ^(r)-c 3 (r)^(,). 



<M>) = o, 



C3(T)^(r)-c 3 (^(i), 



^(^)=^(^)-C3(^(l). 



15 Appendix 3: The first trees 

This appendix gives the Feynman diagrams correspond- 
ing to the first trees for the electron and photon Green 
functions. 



15.1 Electron Green function 

For the electron Green functions, all electron loops are 
oriented anticlockwise and the propagator is oriented from 
right to left, as indicated in <^°(i). Once the Feynman 
diagrams for are known, those for ^(t) are obtained 
by summing all possible insertions of a photon dangling 



bond at each electron propagator — 



and so on for 



(p n (t). Notice that the last two diagrams of <^°(¥) are zero 
by Furry's theorem. However, they are useful to generate 
Feynman diagrams for higher order trees. 



^°(0 
f° (y) 

V°C0 = ■ 
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v ® = ■ ^f ^ 





'AAtCO = ' (J) ' + ' 



+ ■ — Q — • 






. r\ r\ r\ 



15.2 Photon Green function 

For the photon Green functions, all electron loops are ori- 
ented anticlockwise. This is only indicated explicitly for Va m ( i) 
¥>VY). 





VA/i( Y ) = ' — (3 — * 



+ 
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